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We suggest a method for the self-consistent calculations of characteristics of metal films in dielec- 
tric environment. Within a modified Kohn-Sham method and stabilized jellium model, the most 
interesting case of asymmetric metal-dielectric sandwiches is considered, for which dielectric media 
are different from the two sides of the film. As an example, we focus on Na, Al and Pb. We calcu- 
late the spectrum, electron work function, and surface energy of polycrystalline and crystalline films 
placed into passive isolators. We find that a dielectric environment generally leads to the decrease 
of both the electron work function and surface energy. It is revealed that the change of the work 
function is determined only by the average of dielectric constants from both sides of the film. 

PACS numbers: 73.22.Dj, 73.40.Ns, 73.40.Vz, 68.03.Cd 

I. INTRODUCTION 

Thin metal films and flat islands on semiconductor or dielectric substrates can be considered as two-dimensional 
electron systems with properties, which are of interest both from the fundamental point of view and from the per- 
spective of their application in nanoscale electronic devices. 

There are a limited number of experimental works focused on quantum size effects in such systems (for reviews, 
see [H-Q) due to difficulties in sample fabrication, as well as because of lack of suitable experimental methods. One 
of the most important characteristics of metal nanostructures is electron work function. 

There are different methods, which enable one to calculate electron structure of slabs consisting of few monoatomic 
layers. Let us combine them into three groups according to the complexity of computations: I - the Sommerfeld 
electrons in-a-box model (analytical calculations, slabs and wires) [1 Ol — tl 51 ] ; II - self-consistent calculations within 
various versions of jellium model (slabs and wires) [l6l - |20| ; III - ab initio calculations (slabs) [2ll - [24j . The obtained 
results are illustrated in figure 1 for all these three groups. An important ingredient of approaches within group III 
is the monolayer number in the film (see dots in figure 1). For groups I and II, L changes continuously. 

In group I, the Fermi energy (kinetic energy) e-p(L) is counted from the flat bottom of conductivity band, while 
the work function W(L) is counted from the vacuum level. Therefore, their size dependencies are "asymmetric". 
In addition to quantum oscillations, these quantities contain monotonic size contributions, which, at small film 
thicknesses, together show up through inequalities < W(L) < Wq and ep(L) > £p > 0, where Wq and £> correspond 
to the three dimensional (3D) metal (allowing for the energy counting for £>). 

In [H|,[26j], an asymptotic behavior of electron chemical potential for spherical clusters of radius R was determined, 
from which it follows that 

W(R) = W -^<W , (1) 

IX 

ft 2 

where c\ ~ 2.5 eV x o,q for simple metals, ao — r^gy. It is expectable that such a monotonic contribution must appear 
for films also. However, in contrast to the case of group I, self-consistent calculations of groups II and III (see figure 
1), at small film thickness, point out to the suppression of monotonic dependence (having an asymptotic (JTJ) ) by 
corrections of higher orders of smallness. For instance, compensation of terms —jj + jh occurs at L* = and L* 
is large, provided C2^> c\ > 0. 

Experimental results also do not allow to draw unambiguous conclusions on the character of monotonic component 
of W(L): in experiments 0], it is absent (Yb films on Si substrate), while, according to it coincides with the one 

of group I. Note that the comparison of a measured work function for the sandwich consisting of Ag film on Fe(100) 
in 0, [El with calculated results for slabs in vacuum is rather relative. 

Work function as well as film surface energy in dielectric environment have never been calculated before. The 
aim of this work is to compute energy characteristics of metal films in dielectrics. We suggest a method for self- 
consistent calculations of equilibrium profiles of electron concentration, effective potential, energy spectrum, and 
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Figure 1: Illustration of the computation results for groups I, II and III (data for group I are deduced from [lS 
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Figure 2: (a) - Scheme of the film in dielectric environment; (b) and (c) - split semi-infinite metal samples, which have been 
in contact with dielectrics before the splitting. Split parts form a sandwich in figure 2 (a). 



integral characteristics of metal films in dielectrics and dielectric substrates. The developed method is based on a 
stabilized jellium model (27j and local density approximation for exchange-correlation potential (28[, which were used 
by us before [29[ to analyze characteristics of semi-infinite metal with dielectric coating. 

This paper is organized as follows. In Section II, we formulate our model. In Section III, we presents our main 
results and provide a discussion of them. We conclude in Section IV. 



II. MODEL 



Let us consider a metallic film of thickness L at zero temperature. We direct 2— axis perpendicularly to the film 
surface (figure 2 (a), £ ^> L). 

Principal identities for the film can be obtained within a model of a rectangular well for conduction electrons. To 
perform a preliminary analysis, we suppose that the bottom of the potential well is flat and we count energies starting 
from its value. Final expression for the kinetic energies of conduction electrons depends only on energy differences; 
therefore, energies counting in such a way is allowed. 

We study a film of thickness L comparable in magnitude to the Fermi wavelength Af = 27rfcF of an electron 
in 3D metal. The longitudinal sizes of the sample are assumed to be considerably larger than the film thickness 
(L <C L x , L y ), which leads to the pronounced quantization of the transverse component of the electron momentum. 
The three-dimensional Schrodinger equation for a quantum box can be separated into one-dimensional equations. 

The eigenenergies are given by 

kj 

£ik\\ = Ej + — , fc|| = k x + k y , (2) 

where £i is the eigenvalue of the i-th perpendicular state ipi(z) (hereafter the Hartree atomic units are used: H — 
771 = 6 = 1). The eigenvalue e, is the bottom of the i-th subband. For finite and periodic systems in the z-direction 
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Dirichlet and periodic boundary conditions are used, respectively. Therefore, possible allowed electron states k x ,k y , k z 
form a system of parallel planes in the fc— space, k z = fc, (see figure 3). 

Occupation of electron states starts from the point {0, 0, fci} and follows an increase of radius- vector. As a result, 
it turns out that all the occupied states are contained within the area of fc— space, confined between the plane fc 2 = k\ 
and semi-sphere of radius fcF = \/2ef. 

The number of states dZ in each of the circles, formed by the intersection of Fermi semi-sphere with planes k z = kj 
of area S = L x L y , within the interval of wave vectors (fc||,fc|| + dk\\) and taking into account both possible spin 
projections, is 

The maximum value of fcii in each circle numbered by i, is equal to the circle radius fcp(i) = (fcp — ^f) 1 ^ 2 . In order 
to find the number of the occupied states, which coincides with the number of valence electrons N in the film, one 
should integrate dZ over fcn in each circle, and then sum up contributions of all the circles: 

q <-F [■ ' q / «F \ 

* = fx; / *,,*,, = A Uas-x;*? . (3) 

i=l q V i=l / 

Taking into account an expression for electron kinetic energy ^(fc^ + kf), the total kinetic energy of electron 
subsystem equals 

T ° = lt] *n fc n { k i + *i) -it *£ w ( ^ + . (4) 

l=\ Q 1=1 \ / 

where i-p is the number of the last occupied or partially occupied subband. 

In the frame of density- functional theory and stabilized jellium model (SJ), the total energy of metal sample is 
represented by the functional of nonhomogeneous electron concentration n(r): 

£ SJ [n(r)] =T S + E xc + E H + E ps + E M , (5) 

where T s is the (non-interacting) electron kinetic energy, _E XC is the exchange-correlation energy, En is the Hartree 
(electrostatic) energy, £" ps is the pseudopotential (Ashcroft) correction, and Em is the Madelung energy. The sum of 
first three terms in expression ([5]) corresponds to the energy of "ordinary" jellium, Ej. The average energy per valence 
electron in the bulk of metal is £sj, J = Esj,j[n]/N, where N is a total number of free electrons of concentration n, 
defined by valence and atomic density. 

The positive (ionic) charge distribution can be modeled by the step function 



P {z) = n6{L/2-\z\). 



(G) 
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Solving the Kohn-Sham equations 

1 



V ipi (z) + v eS [z, n {z)\ 4>i (z) = Siipt (z) , (7) 



v cS [z, n(z)] = <j>{z) + v xc (z) + (6v) {&ce 9(L/2 - \z\) (8) 
together with the Poisson equation 

V 2 0(z) = -— [n(z)-p(z)], (9) 
e(z) 

with the step function 

C 1; z < -C/2, -L/2 < z < L/2, z > C/2, 
e(z) = I ei; -C/2 <z< -L/2, (10) 
[ e r ; L/2<z< C/2, 

we obtain the single electron wave function and the eigenvalue £j self-consistently. 

We use exchange-correlation potential v xc (z) = d[n(z)e xc (z)]/dn(z) in LDA(28|. The term (5v)f acc in ©, which 
makes it possible to distinguish different crystal faces, represents the difference between the potential of the ionic 
lattice and the electrostatic potential of the positively charged background averaged over the Wigner-Seitz cell: 



(£«)face = (Sv)wS ~ + i ( fo )wS = -n — 



V 3 b ) dn 

where d is the distance between the atomic planes parallel to the surface. The term (5v)yfs describes a polycrystalline 
sample [13] ■ In equation (|10l) e\ and e r are dielectric constants of isolators from the left and right side of the film, 
respectively. 

The electron density profile n(z) is expressed through the wave functions ipi(z) 
Values of %y and £f are determined by the solution of the equation 

i P 

i-pSY — irLh + e,-; i — 1,2, ... , ip; Si < £p, (12) 

i=i 

which follows from the normalization condition ([3]) and definition of the Fermi energy. 

In nanofilms, Friedel oscillations are significant throughout the sample. Therefore, energies are counted from the 
vacuum level, which is the energy of the electron in rest in the area \z\ ^> C/2. For bound states, energies are negative, 
including £p. 

We use iterative procedure (see Appendix A) allowing us to solve self-consistently the system of equations ([ZJ) , ([9]) , 
(JTTj) and to find optimal profiles n(z), 4>(z), as well as spectrum of one-particle energies. As a result, work function 
is defined in the form 

W = -s F . (13) 

Let us consider a scheme for the surface energy determination (see figure 2 (b) and (c)) for the film of thickness L 
in a dielectric environment. 

First, we take a semi-infinite metal (Meoo) covered by a dielectric (e r ). Let us denote the energy of such a sample 
as i?{Meoo|er}- We now split the sample and move the parts, as shown in figure 2 (b). As a result, two new surfaces 
of the same area S are formed, which are in a contact with the vacuum (e = 1). We denote the energies of these two 
parts as ^{Meooll} and £?{l|Mc£/ 2 | e r}, while the irreversible work A, which is needed to form them, as 

£{Mcoo|l} + £{l|Me L/2 |e r } - E{M eoc \e r }. (14) 

Let us stress that, as a result of these manipulations, the "fabricated" sandwich represents a film of thickness L/2 > d 
on the dielectric substrate in vacuum (air). 
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Similar manipulations with another sample (figure 2 (c)) require a work 

E{l\M Coo } + £{ei|Me i/2 |l} - £{ ei |M eoo }. (15) 
It is convenient to represent the total energy as a sum of bulk and surface contributions 

E = E h +E S . 

Then, bulk components E h do compensate in the expressions (fT4|) and (fT5j) . In each case considered above, the specific 
surface energy 7 equals A/2S. 

The work needed to "create" a film on a dielectric is 



A{e|Me i/2 |l} = i [^ s {e|Me i/2 |l} + E s {l\M eoo } - E s {e\M eoo } 



(16) 



Now, we join two sandwiches by their free surfaces. We then obtain a film shown in figure 2 (a). The work to create 
it can be represented as the energy of adhesion of such two pieces with the minus sign 



A{ei|Mei|e r } = i [^ 8 {ei|Mei|e r } - £ s {ei|Me i/2 |l} - £ s {l|Me i/2 |e r } 



(17) 



Electron density profiles and potentials for each of the contributions in the expressions (fl6|) and (fl7|) are different, so 
that they must be calculated self-consistently and separately. 

As similar to the definition for the semi- infinite metal [2!| l3l| . E s for the film is determined by the difference 
between the total film energy ([5]) and the energy of homogeneous metal (stabilized jellium) of the same volume: 



L/2 



-L/2 



E s {e l \Me L \e r }^Esj(L)^SLnesj^2s!^fj + (Sv)! llcc J dz[n(z) 
By using equation (jj) , quantum-mechanical definition of an energy 



(18) 



as well as the definition given by equation (|18p . we obtain an expression for the first component of 7j: 

1 



7s 




dzil)i(z)W ipi(z) 



-Lnt« 



where t s = 3fcp/10 is the kinetic energy per 1 electron for bulk. The remaining components are 



dz n(z)e xc [n(z)] - -Lne xc (n); 



(19) 



(20) 



7H = i J dzcf>(z)[n(z)-p(z)}. (21) 
— 00 

For asymmetric sandwiches, {ei|MeL|e r }, due to the formal division on the doubled area, the surface energy is cal- 
culated "in average" . This is the consequence of the definition of 7 through the integral of tangential component of 
pressure tensor over z from — 00 to +00. The pressure tensor contains the nonelectrostatic part and the Maxwell 
stress tensor (cf e.g. [Hj]). 
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III. RESULTS AND DISCUSSION 



We perform calculations for both polycrystalline and crystalline films made of Na, Al and Pb, with electron 
concentration n = 3/4irr^ with corresponding electron parameter r s = 3.99, 2.07 and 2.30 ao. 

The minimal thickness of "crystalline" sandwiches should be not less than 2d. It must be equal to 4d for {ei|Me|e r }, 
only in the case Eq. (fTT]) is used, d is comparable to ^Af (Af = 13.06, 6.78 and 7.53 ao for Na, Al and Pb, 
respectively). 

In view of the multimolecular thicknesses of dielectric coatings on the metal film surfaces and rapid fall of the electron 
distribution outside of a film (approximately at a distance of 10 - 15 ao), we formally neglected the effect of a thickness 
of the coatings, whose minimum thicknesses must be much greater than that of a monatomic (molecular) layer of a 
dielectric. For comparison, we recall that the diameter of atom Si equals 5 ao. If the thicknesses of the dielectrics 
on the surfaces of the metal film are foreseen to be less than the free path of an electron (tens of angstrom) injected 
into these dielectrics, then the calculated work function is qualitatively appropriate for the vacuum-dielectric-metal 
slab-dielectric- vacuum. 

For symmetric sandwich the effect of a dielectric coating on the surfaces is reduced to the "elongation" of the 
electron distribution tail and the effective potential beyond the surface of a metal (polycrystalline films {1|A1|1} and 
{3|A1|3} on Fig. 4). The calculations were performed for e = 1,..., 12. Inside the film one can see the Friedel 
oscillations of electron density with peaks near geometrical boundaries. The period of oscillations is close to ^Af and 
only weakly depends on the presence of dielectric coatings. The situation is similar for Na and Pb films. 

At the boundaries between the metal film and the coatings, there are jumps in the derivative of the electrostatic 
potential 4>'(z), which disappear, provided the dielectric constants of the coatings are equal to 1. These jumps are 
due to the boundary conditions (|A2[) at z = ±L/2. The jumps are also reflected on v e s(z) profile, since <p(z) is one 
of its components. In addition, at the borders, there are another jumps of not only the derivative v' eS (z) 7 but also of 
v e ff(z) profile itself for any values of e, including e = 1. Such jumps have another origin compared to the first ones. 
This fact is linked to some features of the model [27j, namely to the presence of the effective potential component 
(<5i>) face 6(L/2 — \z\). These nonphysical jumps should not be taken into account in the estimation of the effective force 

It is seen from figure 4 that force orientations are opposite at both sides of the film, so that the film in whole must 
be stressed. The existence of the force should lead to the increase of spacings between some lattice planes d, while 
spacings between other planes must become narrower. 

The depth of the potential well \v c g\, in which the electrons are located in metal film, decreases "in average" with 
increasing e and, as a result, the electron work function also decreases (see figure 6). 

Film spectra {l|Ali|l} are presented in figure 5. For comparison, in the same figure, we also provide the results 
obtained within the electrons-in-a-box model with the well depth Uo = — (Wo + £f) < 0. 

It is seen from figure 5 that the dependence of the eigenstate energies on the film thickness, within the SJ model, is 
oscillating and decreasing. For subbands with large numbers i = 10, 11, there are gaps due to the algorithm instability 
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Figure 5: Results of calculation for the energy spectrum (subbands) and Fermi energy ef{L) of the film {1|A1|1} by the 
self-consisting method (solid lines) and in rectangular-box model (dashed lines). 

in the vicinity of the vacuum level. Within the rectangular-box model, this dependence is only decreasing. Due to 
smoother edges of the self-consistent well, it contains more subbands compared to the model of a rectangular box. 
Difference in subbands numbers significantly affects calculated dielectric function and optical conductivity of the 



Within the rectangular-box model, in contrast to the SJ model, e-p(L) is always located above one for 3D metal. 
Amplitudes of oscillations decrease as L increases. Within both models, maximum Fermi energies (minimum work 
functions (1131) ) correspond to the points, in which curves of eigenenergies intersect Fermi energies. Within the SJ 
model, in contrast to the rectangular-box model, minimum Fermi energies correspond to the points, in which Fermi 
energy is located between two nearest eigenenergies (magic film thicknesses similar to magic numbers in clusters). 

Asymmetric sandwiches {ei|Me|e r } and {l|Me|e}, which contacts the air or vacuum, are of particular interest from 
the viewpoint of experimental investigation due to the perspective of their use in technological applications. For 
instance, these are films on Si (e = 12) or Gc (e = 10) substrates with different concentration of impurities and 
crystallographic orientations (see, e.g. |5|). 

Let us consider both electron density and potential profiles for the polycrystalline film {1|A1|5}. Presence of a 
dielectric at the right side of the film leads to the asymmetry of electron distribution (see the insets in figure 4), so 
that there appears a hump in both the electrostatic and effective potential at the left side above the vacuum level. 
This should result, for example, in the anisotropy of a field emission along the z— axis. It is worth mentioning that 
bottoms of wells for sandwiches {1|A1|5} and {3|A1|3} are essentially the same, some difference appears only in "tails" 
of potential profiles. 

It is of interest to compare heights of humps at L = 10, 12, 13.5 and 20, 22, 23.5 ao- These thicknesses correspond 
to the minimum, maximum, minimum of the dependence W(L) for {1|A1|5}. It turns out that, with the increase 
of L, the hump height weakly oscillates and decays similarly to the work function, but maxima of the hump height 
corresponds to minima of the W(L). For the values of L, as given above, these heights are 0.176, 0.148, 0.170 and 
0.158, 0.139, 0.156 eV, respectively. 

In order to analyze such a behavior of potential profiles, it is necessary to go beyond the isotropic model based 
on a defined ([5]) distribution of homogeneous positively charged background, i.e. one has to take into account not 
only the reaction of the electron subsystem, but also the reaction of the ion subsystem to the presence of a dielectric. 
Spacings between the lattice planes are determined by the balance of forces from the right and left sides for each 
plane. A simplest realization of this idea is to disregard variations of spacings between the lattice planes and to vary 
the profile of the ion jellium distribution ©. We found that such a procedure leads to a significant deformation of 
the well bottom, but does not result to considerable changes of both the spectrum and hump height. 

Figures 6 show results of our calculations of both the electron work function and surface energy for crystalline 
sandwiches using expression (fTE)) . Horizontal lines correspond to semiinfinite samples. In contrast to the surface 
energy, size dependences W(L) have deep and pronounced minima. It is easier to analyze them using a simple model 
Amplitudes of largest work function "oscillations" are smaller than 0.5 eV. By considering dependencies for 
different metals, it is easy to see that all the differences are due to values of r s . For the Al, which has the smallest r s , 
work function oscillations are maximum, while the period is minimum. Positions of both maxima and minima depend 
weakly on e of a dielectric and slightly shift towards smaller L with the increase of e. 
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Figure 6: Work function and energy per unit of area for crystalline sandwiches {ei|Me|e r } and semi-infinite metal covered by a 
dielectric {Meoo|e} (Me = Na, Al, Pb). 



In contrast to the work function, surface energy oscillations can be approximated by analytical dependences 

sin (2kpL + if) 



E s {e\Me L \e\ = E s {M Coo \e} + A- 



L 



with parameters A and if. Maxima of function W{L), "f(L) correspond to "magic" film thicknesses, which are defined 
by maximum occupation of a given subband. 

The unexpectable result of self-consistent calculations is a coincidence of dependencies W(L) for sandwiches 
{l|Me|12} and {6.5|Me|6.5}. Computations for {l|Me|5} and {3|Me|3} give a similar result. This means that the 
electron work function for asymmetric sandwiches {ei|Me|e r } coincides with high accuracy with the work function for 
symmetric sandwiches {(e)|Me|(e)} with the averaged value (e) = + e r ). 

Work function has both the bulk and surface contributions. Because bulk metal contributions W(L) for sand- 
wiches {l|Mei|12} and {6.5|MeL|6.5} are the same by definition, also the same are contributions of dipole surface 
barriers. We here imply the total contribution of both sides of a sandwich, since the work function is an "isotropic" 
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characteristics 33] . A coincidence of work functions is most likely a geometric effect. This feature will be addressed 
elsewhere. 

For surface energies, such a coincidence does not exist. It is not difficult to perform calculations according to 
formulas (|I6I) and flTTl) . if 7 are known. 

IV. CONCLUSIONS 

We proposed a method for the self-consistent calculations of spectra, electron work function, and surface energy of 
metal films placed into passive dielectrics. As typical examples, we considered Na, Al, and Pb films. 

The effective force acting on the film from the outside is due to the inhomogeneous electron distribution. This force 
should lead to film stressing in a transverse direction. The effect of the stressing generally becomes more significant 
with the increase of the film thickness. 

In contrast to the surface energy, size dependencies of work function have deep and strongly pronounced minima. 
The smaller r s the more difficult the problem of numerical analysis of size dependencies in the vicinities of these 
minima. 

With the increase of film thickness up to few Af, size variations of both the work function and surface energy 
occur near their average values (for symmetric sandwiches, these values correspond to 3D metals and do not contain 
significant monotonous size contributions). Dielectric environment generally leads to the decrease of electron work 
function and surface energy. 

We also considered asymmetric metal-dielectric sandwiches characterized by different dielectrics at both sides of 
the film. One of the examples of such systems is a film on the dielectric substrate. We found that the presence of a 
dielectric from one side of the film leads to such a "deformation" of electron distribution that there appears a "hump" 
above the vacuum level both in the electrostatic and effective potentials. The asymmetry of potential profile of this 
kind should lead to an anisotropy of the field emission. In addition to size dependencies, the shift of the work function 
is generally determined by the average dielectric constants of environments. 

Obtained wave electron functions for metal-dielectric sandwiches allow us to calculate the matrix elements of optical 
transitions, conductivity tensor, and coefficient of optical absorption fl4|. 
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Appendix A: Self-consistency procedure 



The initial approximation n(z) is chosen for solving the Kohn-Sham equations in the form of a one-parametric trial 
function n^{z) = nf(z), where 

( _l e (*-L/2)/A + l e (*+L/2)/A j z < _ L f 2i 

f(z) =\l- leM- L ' 2 )/ x - i e -(*+£/2)/A \ z \ < L/2, 
[ _| e -(»+£/2)/A + i e -(«-L/a)/A j z > L/2 . 

A is the variational parameter, which is found through the minimization of surface energy. Solution by a direct 
variational method is an independent problem, which is not addressed in this paper (for simple metals A is closed to 
1 ao). As a result of integration of equation (jH]), within the initial approximation, we obtain <j>(°\z) = — 47mA 2 /(z). 
Each wave function ip(z) is constructed as 



\ fright (z), Z > Zq, 



under the condition of continuity of functions ipi e it(zo) = Vvight^o), as well as of their derivatives ip[ e{t (zo) = ip' light (zo) . 
zq is an arbitrary point in the interval z € [—L/2; +L/2], while ^left^) and bright (z) are functions, which are found 
by a numerical solution of Eq. ([7]) by the Numerov's method from z = z_ to z = z$ and from z = z + to z — zq, 
respectively. It is sufficient to take values z T — ^f(L + 20) ao. In these points, the potential profile v e g(z) is cut off. 
The boundary conditions (J7]) here are determined by the behavior of the wave function ifj under the barrier from the 

left (e z V^^) and right (e^V^^) sides from the slab (|z| > | | ) , respectively. Boundary conditions provide wave 
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function, as well as its derivative at z — z T . This peculiarity of our computations is due to the fact that errors of 
the numerical method for the wave function bright (z) and ip\ c ft(z) near the right and left boundaries of the interval 
grow, since the round-off errors also increase and lead to the instability of the algorithm under the motion towards 
the exponential damping. 

In order to solve the system of equations ([7]) , (O and (fTTj) self-consistently, with relatively small number of iteration 
steps, the Poisson equation ([9]) should be modified, in particular, by introducing a perturbation (30j |. 
Equation @ is solved by the Lagrange method in the form 



,2^,0-1) _ 



47T 



rft' — p 



(Al) 



with the boundary conditions 



^t(z) = ^{z), e^£ t '(z) = <f$\z); z = -L/2, 
= 4>£t(z), 4$\z) = e r 0«? t '(z); z = L/2, 
0; z — =foo. 



(A2) 



The term q 2 <p was introduced as a small perturbation; cf> ut(z) an d ^in(^) are potentials outside and inside the film, 
respectively. In equation (|A1|) . at each step of the iteration j — 1, 2, 3, electrostatic potential profile depends not 
only on the electronic concentration profile, but also on its own profile at the previous iteration. It is convenient to 
take q equal to electron wave number at the Fermi sphere k-p = (S^n) 1 / 3 of homogeneous electron liquid. 
The solution of equation (|A1[) for C — >• oo has the simple form 



■ U \z) = I 



j ^-fidz' + A 1 )e^+[ - J ^fidz' + B 1 ) e-«*, z < -L/2 



2q 



J <^-f 2 dz' + A 2 e'* + - / ^f 2 dz' + B 2 e-«', |*| < L/2, 

-L/2 / V - L / 2 / 

OO _ / \ / OO / \ 

- I ^-fadz' +Aa)e*>+(J ^-fsdz' + B 3 e~*; z > L/2, 



(A3) 



where f m (z') — — An[n(z') — p(z')]D m — q 2 ^" 1 ^ (z ! ) and D m = ej~ , 1, e^ 1 for m = 1,2,3, respectively. The choice of 
values B\ = and A3 = immediately follows from the condition of finiteness of potentials far away from the film. 
Values of other coefficients A and B are found from the solution of the system of equations (jA2l) : 



Ai = 



-L/2 

2A 2 1-ei f e<^'+ L ) 



-L/2 



1 + ex 1 + ex 



2g 



/id*' - 



2q 



-fidz', 



B* = 



2B, 1 



L/2 



1 + £r 1 + £ r J q 

-L/2 



, 1 - e r /" e 



,-q(z -L) /■ ijz 

1 + e r 7 2g 7 2g 

L/2 L/2 



Let's introduce notation 



J{±) = Yo 



-L/2 L/2 L/2 

Fl2q(lTer) / <fe'e"7i + Y 2 (l ± q)(l + e r ) / dz'e- 92 ^ + Y 3 {1 ± q)(l - e r ) J dz'e^'h 



-L/2 



-L/2 
00 

r 4 2e r (l± ei ) /" dz'e-^'h 

L/2 



, (A4) 



where F = {2g[(l - ei)(l - e r )e-« L - (1 + ei)(l + e^]}" 1 . Then A 2 = J (+) for Y 1)3 = 1, Y 2A = e« L and B 2 = J ( _ } 
for F 2 ,4 = 1, Y x = e« L = F 3 _1 . 
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In the case of the symmetric sandwich e\ = e r the accurateness of calculations is verified by examination the 
stationarity conditions n'(z) = and <f>^ {z) = in the center of the slab (z = 0). 
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